Is A^-gap-only superconductivity possible in Mgi_ x Al x B 2 and Mg(B\- y C^)2 alloys ? 
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Using density-functional-based method, we study the k-resolved u- and 7r-band holes in 
Mgi- x Al x P>2 and Mg(Bi- y C v )2 alloys. We find that the calculated profiles of the loss of a- and 
7r-band holes in these two systems as a function of impurity concentration are in qualitative agree- 
ment with experiments, as expected. We also describe its implications vis-a-vis superconductivity 
in Mg 1 _ x Al x B 2 and Mg{B 1 _ y C y ) 2 . 



The surprising discovery of superconductivity at 40 K 
in a simple metal MgB 2 [l| a few years ago contin- 
ues to test the details of the microscopic picture offered 
by Bardeen-Cooper-Schriefer theory of superconductiv- 
ity and its extensions thereof. Generally, it is agreed 
that in MqB 2 (i) there are two superconducting gaps 
H H II H U B, A CT and A T of 7 meV and 2 meV, re- 
spectively, (ii) the two gap s arise from the holes in the B 
Per and Pit bands [8l lit lift llll ll2L ll-'l Il4|. (iii) the holes 
in the B p„ and p v bands coup le strongly to the in-plane 
E 2g phonon mode [lj, LLij, |lg, LLD] , and (iv) the coupling 
is particularly strong along T-A direction and around M 
point in the hexagonal Brillouin zone 

To gain further insight into the role played by the holes 
of the B p„ and p^ bands in deciding the superconduct- 
ing properties as well as the relative strengths of the two 
gaps A a and A v in MgB 2 , recently, there has been a 
concerted effort to examine the various superconducting 
properties of doped MgB 2 . In particular, the Al dop- 
ing at the Mq site and the C doping at the B site in 
MgB 2 [13, E3, EH E3, (S3, E3, ISH [2a E3, E3, IS^I are 
expected to fill up the holes in the B p a and p v bands 
which, in turn, would affect the two gaps A a and A v . 
Understanding the changes in the superconducting prop- 
erties of Mgi_ x Al x B 2 as a function of Al concentration 
x and of Mg(Bi^ y C y ) 2 as a function of C concentration 
y would provide a more detailed information about the 
electron-phonon and the electron-electron interactions in 
these alloys. 

Experimentally, in Mg\- x Al x B 2 the superconducting 
transition temperature T c as well as the gap A „ ar e found 
to decrease with increasing Al concentration l2ll I22L 
151 The change in the gap A^ is relatively small 

up to x 0.7 at which the superconductivity vanishes in 
Mch-f ALBo. However^ in Mg(B 1 - y C v ) 2 the T c filial 
EllHl^EaESEllilElliDl and the A^ decrease rapidly 
as a function of C concentration with superconductivity 
vanishing at around y 0.15. The change in A^, once 
again, is found to be minimal in Mg(B\_ y Cy) 2 [2I 151 



Most of the previous theoretical studies jd.l 144 Idol l4(l 
H El E3 IH3 of Mg^ x Al x B2_and_ Mp(Bi-«C y ) 2 
used virtual-crystal app roximation |4.T . I47L l4?4 l50L l5l| or 

supercell approach [46], having limited predictive capa- 
bility, to describe the chemical disorder in the Mg sub- 



lattice and the B sub-lattice respectively. An accurate 
and reliable description of chemical disorder, provided 
by coherent-potential approximation, has been used pre- 
viously in the study of Mg\- x Al x B 2 and Mg{B\_ y Cy) 2 
[44L 14.4 l4fl | but these studies provided either k- integrated 
information or were inconclusive. 

Since many of the superconducting properties of MgB 2 
and its alloys Mg\_ x Al x B 2 and Mg(B\_ y Cy) 2 are cru- 
cially dependent on their normal state electronic struc- 
ture 18, 15, 17, IHl IHl IHl ] . especially along T-A direction 
and M point of the hexagonal Brillouin zone (BZ), we 
study the changes in Mg\- x Al x B 2 and Mg{B\- y Cy) 2 al- 
loys as a function of Al and C concentrations respectively, 
using k-resolved Bloch spectral function A(k, E) (HllHg 
and B and C p x ( y ) and p z densities of states. It is ex- 
pected that with increasing Al and C concentrations the 
holes in the respective alloys would get gradually filled 
up. However, not surprisingly, the way the holes are filled 
up in Mg\- x Al x B 2 and Mg(B\- y Cy) 2 alloys make the 
two systems quite different from each other as detailed 
below. 

In this Letter we show that in Mg\_ x Al x B 2 (i) the ri- 
band holes along T-A get annihilated by x 0.8 with the 
holes at T vanishing first at a; 0.5, and (ii) the 7r holes 
at M vanish by x 0.7. However, in Mg(Bi^ y C y ) 2 (i) 
the er-band holes along T-A deplete rapidly and vanish 
beyond y ~ 0.25 with the holes at A vanishing first, 
and (ii) the 7r-band holes at M as well as a-band holes 
at A reduce substantially beyond y ~ 0.1. Before we 
discuss our results, we briefly describe the details of our 
calculations. 

The normal metal electronic structure of the disor- 
dered alloys Mg\- x Al x B 2 and Mg(Bi- y C y ) 2 are calcu- 
lated using the Korringa-Kohn-Rostoker (KKR) Green's 
function method formulated in the atomic sphere approx- 
imation (ASA) (HIIH^. which has been corrected by the 
use of both the muffin-tin correction for the Madelung 
energy, needed for obtaining an accurate description of 
ground state properties in the ASA and the multi- 
pole moment correction to the Madelung potential and 
energy which significantly improves the accuracy by tak- 
ing into consideration the non-spherical part of polar- 
ization effects 60]. Chemical disorder was taken into 
account by means of coherent-potential approximation 
(CPA) [fil| . The exchange and correlation were included 
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within the local density approximation (LDA) using the 
Perdew-Wang parameterization of the many-body cal- 
culations of Ceperley and Alder jfi^]. During the self- 
consistent procedure the reciprocal space integrals were 
calculated by means of 2299 k-points in the irreducible 
part of the hexagonal BZ, while the energy integrals 
were evaluated on a semicircular contour in the com- 
plex energy plane using 24 energy points distributed in 
such a way that the sampling near the Fermi energy was 
increased. When calculating the density of states the 
number of k-points were increased to 6075. The atomic 
sphere radii of Mg and B were kept as 1.294 and 0.747, 
respectively of the Wigner-Seitz radius. The sphere radii 
of the substituted Al and C were kept the same as that 
of Mg and B respectively. The overlap volume resulting 
from the blow up of the muffin-tin spheres was approx- 
imately 10%, which is reasonable within the accuracy 
of the approximation The charge self-consistency 

iterations were accelerated using the Broyden's mixing 
scheme j^J], the CPA self consistency was accelerated 
using the prescription of Abrikosov et al |§3 ■ The calcu- 
lations for Mg\^ x Al x Bi for < x < 1 and Mg{Bi_ y C y ) 2 
for < y < 0.3 were carried out at experimentally 
obtained lattice parameters of these alloys , assuming 
a rigid underlying lattice. The effects of local lattice- 
relaxation as well as any possible short-range ordering 
effects are not considered. Our results for Mg\- x Al x B 2 
and Mg(B\-yCy)2 are analyzed in terms of Bloch spec- 
tral density, j4(k, E) j^| and the B and C p densities of 
states. 

In general, electron doping of MgB 2 by substituting Al 
at the Mg sub-lattice or C at the B sub-lattice can have 
three important effects, (i) the outward movement of the 
Fermi energy due to addition of electrons, (ii) disorder- 
induced broadening of peaks in A(k, E), and (iii) redis- 
tribution of states due to hybridization. If the disorder 
and hybridization-related effects are negligible then elec- 
tron doping can be well described by a rigid-band-like 
picture, where Ep is shifted up using MgB 2 densities of 
states until the added electrons are accommodated. In 
such a case, the a- and 7r-band holes would gradually get 
filled up as a function of electron doping. Thus, in the 
rigid-band picture the a- and 7r-band-hole-dependent su- 
perconducting properties of MgB 2 alloys would be sim- 
ilar as a function of electron doping. However, if the 
disorder and the hybridization lead to significant modi- 
fications of the a- and 7r-band holes then the electronic 
properties of both the normal and the superconducting 
states of MgB 2 alloys would depend on the details of 
electron doping rather than just the electron count. Ex- 
perimentally, as outlined in earlier paragraphs, the a- 
and 7r-band-hole-dependent properties of the two alloys, 
Mgi- x Al x B 2 and Mg{B\- y C y ) 2 , are found to differ sig- 
nificantly from each other as a function of Al and C con- 
centrations, underscoring the importance of disorder and 
hybridization effects in these alloys, in particular on a- 



and 7r-band holes. Therefore, in the following, we pri- 
marily focus on the changes in the states close to Ep, 
which happen to be the p a states at T and A points and 
p,r state at M point of the hexagonal BZ. 




Figure 1: The Bloch spectral function A(k,E) of 
Mgi- x Al x B2 alloys calculated at F (top panel), A (middle 
panel) and M (bottom panel) points of the hexagonal Bril- 
louin zone using the KKR-ASA-CPA method for x equal to 
(black), 0.1 (red), 0.3 (green), 0.5 (blue), 0.7 (yellow), 0.9 
(brown) and 1 (cyan), respectively. The vertical line through 
the energy zero represents the Fermi energy. 

The most important result of the present work for 
Mgi^ x Al x B 2 alloys is shown in FigQ] where we show 
the Bloch spectral function A(k, E) as a function of x for 
< x < 1 at r, A and M points of the hexagonal BZ in 
a limited energy window (±0.1 Ry) around Fermi energy 
Ep. Along T-A, the doubly degenerate B p a states in 
Mg\_ x Al x B 2 alloys move closer to Ep with increasing 
Al concentration from x = 0.1 to x — 0.70, indicating a 
decrease in er-band holes and subsequently, weakening of 
electron-phonon coupling along this direction. We also 
find that for x ~ 0.5, the peak in ^4(k, E) is at Ep indi- 
cating some topological changes in the Fermi surface of 
the alloy. By x — 0.90, the B p a states at A are well 
inside Ep , resulting in the annihilation of all the p a hole 
states along T-A. Thus, we should not expect er-band- 
hole-based superconductivity beyond x = 0.80. We also 
find that the Mg-derived antibonding s-state (Yq) moves 
closer to the p a state for intermediate concentrations but 
separates out for x = 1.0. 

The B p,r hole states at M in Mg\- x Al x B 2 alloys get 
gradually annihilated with increasing Al concentration 
from x = 0.1 to x = 0.50. We find that the A(k, E) 
peak corresponding to p^ states crosses Ep at a lower x 
than the p a states at A point. As the contribution to the 
electron-phonon coupling responsible for 7r-band super- 
conductivity comes from M and other similar points in 
the BZ, a A^-gap-only superconductivity seems unlikely 
in Mg\^ x Al x B 2 alloys. However, it does not preclude 
a crossover from a a to 7r dominated superconductivity 
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2: The Bloch spectral function A(k, E) of 
yCy)2 alloys calculated at Y (top panel), A (middle 
panel) and M (bottom panel) points of the hexagonal Bril- 
louin zone using the KKR-ASA-CPA method for y equal to 
(black), 0.05 (red), 0.1 (green), 0.15 (blue), 0.2 (yellow), 
0.25 (brown) and 0.3 (cyan), respectively. The vertical line 
through the energy zero represents the Fermi energy. 



In Fig. El we show the Bloch spectral function A(k, E) 
of Mg{B\- v Cy)2 alloys as a function of y for < y < 
0.3 at T, A and M points around Fermi energy. We find 
that with increasing C concentration from y = 0.05 to 
y = 0.2 the p a states along T-A move towards Ep as 
well as get redistributed on the energy scale due to dis- 
order and hybridization. In fact, in the range y = 0.1 
to y = 0.15, we find substantial reduction in p a states 
around A point. Since the electron-phonon coupling is 
stronger near A point, the loss of p a states should dra- 
matically affect the cr-band superconductivity. Similar to 
the reduction of p a states along T-A, the p^ hole states at 
M point in Mg{B\- y C y )2 alloys are drastically reduced 
by y — 0.1. Here, it must be noted that the reduction in 
a- and 7r-band holes in Mg\- x Al x B2 with increasing Al 
concentration is mostly due to the upward movement of 
the Fermi energy. In contrast, disorder and hybridization 
play a dominant role in reducing the a- and 7r-band holes 
in Mg(Bi- y Cy)2 alloys. 

The changes in the B p a and p^ hole states in 
Mgi_ x Al x B2 as a function of Al concentration and in 
Mg(Bi- y Cy)2 as a function of C concentration can be 
seen more clearly in their respective partial densities of 
states as shown in Fig. In Mg\- x Al x Bi alloys the 
B Px(y) hole states decrease steadily as a function of x, 
vanishing for x ~ 0.8 as shown in Fig. 03 Thus, as stated 
earlier in the context of A(k,E), we should not expect 
cr-band-hole-based superconductivity beyond x = 0.80 
Mgi- x Al x B2- However, the decrease in B p z hole states 
with increasing x is not as much. Even for x — 1.0, there 
are B p z hole states available in Mg%- x Al x B2 as shown 
in Fig|21 but it does not sustain A„.-gap-only supercon- 




Figure 3: The atom (B and C) and orbital (p^^and p z ) 
projected density of states (DOS) in Mg 1 ^ x Al x B 2 (left pan- 
els) and Mg(Bi- y C'y)2 (middle and right panels) alloys cal- 
culated using the KKR-ASA-CPA method for x(y) equal to 
0(0) (black), 0.1(0.05) (red), 0.3(0.1) (green), 0.5(0.15) (blue), 
0.7(0.2) (yellow), 0.9(0.25) (brown) and 1(0.3) (cyan), respec- 
tively. The vertical line through the energy zero represents the 
Fermi energy. 



ductivity. 

For Mg(B\_yCy)2 alloys, our results show that p x ( y ) 
and p z densities of states of both B and C behave sim- 
ilarly with increasing C concentration, as can be seen 
from Fig. By y ~ 0.25, most of the p x ( y ) hole states of 
B and C have moved inside Ep. As pointed out above, 
in Mg(B\_yCy)2 alloys with the addition of C those cr- 
band holes are lost first which couple to the phonons 
more strongly which may lead to a loss of superconduc- 
tivity in spite of having some p x ( y ) hole states of B and 
C. The changes in the p z hole states are minimal in 
Mg{B\- y Cy)2, with a slight increase in the hole states 
for higher C concentrations. 

We like to point out that the approximations used in 
the present work, namely the use of spherically sym- 
metric potential to describe the effective one-electron in- 
teraction and the overlapping atomic spheres to fill the 
space, are expected to affect some numerical aspects of 
the present work. The use of optimized lattice constants 
is expected to similarly affect the present work, especially 
in Mg(Bi- y Cy)2 alloys for y > 0.2 because for these 
concentrations we have used the lattice constants corre- 
sponding toy = 0.15. However, our results for MgB2 and 
AIB2 are consistent with the more accurate, full-potential 
results, and since the coherent-potential approximation is 
known to describe the effects of disorder reliably, we ex- 
pect our intermediate-concentration-results to be robust 
and qualitatively correct. 

In conclusion, we have studied the k-resolved a- and 
7r-band holes in Mg\_ x Al x B2 and Mg{B\- y Cy)2 alloys. 
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The calculated profiles of the loss of a- and 7r-band holes 
in these two systems as a function of impurity concentra- 
tion are expected to be in qualitative agreement with 
the experiments. We have also described its implica- 
tions vis-a-vis superconductivity in Mg\— X M X B<2. and 
Mg{B]_- y Cy)i alloys. 
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with Dr. Igor Mazin on the electronic structure proper- 
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